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Abstract. The CeMIns heavy fermion compounds have attracted enormous interest 
since their discovery six years ago. These materials exhibit a rich spectrum of unusual 
correlated electron behavior, and may be an ideal model for the high temperature 
superconductors. As many of these systems are either antiferromagnets, or lie close to 
an antiferromagnetic phase boundary, it is crucial to understand the behavior of the 
dynamic and static magnetism. Since neutron scattering is difficult in these materials, 
often the primary source of information about the magnetic fluctuations is Nuclear 
Magnetic Resonance (NMR). Therefore, it is crucial to have a detailed understanding 
of how the nuclear moments interact with conduction electrons and the local moments 
present in these systems. Here we present a detailed analysis of the hyperfine coupling 
based on anisotropic hyperfine coupling tensors between nuclear moments and local 
moments. Because the couplings are symmetric with respect to bond axes rather than 
crystal lattice directions, the nuclear sites can experience non- vanishing hyperfine fields 
even in high symmetry sites. 
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1. Introduction to NMR in the CeMIns Materials 

The CeMIns materials (also known as the 115's), with M= Co, Rh or Ir, exhibit a 
rich variety of phenomena as a result of strong electron-electron correlation . These 
phenomena include unconventional d-wave superconductivity, local moment magnetism, 
coexistence of magnetism and superconductivity, and a possible Fulde-Ferrell-Larkin- 
Ovchinnikov superconducting phase OIHl- The superconducting transition temperature, 
Tc, of CeColns is the highest of the Ce-based heavy fermion materials. Recently, an 
isostructural class of plutonium based 115's was synthesized, the PuMGas series, with 
M= Co or Rh, that are superconductors with T^s of 18. 5K, and 8.5K, respectively |3]. 
Since PuCoGas is a d-wave superconductor like CeCoIus, it has been argued that its ten- 
fold increase in is due to an increase in the size of the relevant magnetic interaction 
energy, J, and thus the PuMGas materials bridge the gap between two distinct sets of 
unconventional superconductors: the heavy fermion superconductors, with Tc~ 1 — 2K, 
and the high temperature superconductors, with Tc~ lOOK [S]. Consequently, one can 
reasonably argue that by investigating the role of magnetism and superconductivity in 
the CeMIns materials, one can gain considerable insight in the the role of magnetism and 
superconductivity in the high-Tc materials and perhaps shed new light on the details 
of the pairing mechanism. These heavy fermion superconductors allow us to probe 
the physics over a much broader temperature range compared to the characteristic spin 
fluctuation energy, J (where J typically is on the order of 1500K for the cuprates and lOK 
for the heavy fermions). Furthermore, the CeMIns materials are more easily handled in 
the laboratory, can be tuned with pressure rather than doping, and are extremely clean 
with small residual resistivities and large scattering lengths. 

The magnetism in the CeMIns system is clearly relevant at some level for the 
existence of the unconventional superconductivity and the non-Fermi liquid behaviors 
jni El El E]- The magnetism can be investigated experimentally by measuring the 
dynamical magnetic susceptibility experimentally, either with neutron scattering or 
NMR. Neutron scattering experiments often are limited because the ^^^In strongly 
absorbs neutrons (cross section ~ 200 barns), which severely limits the resolution. On 
the other hand, NMR can provide important insight into the behavior of the dynamical 
susceptibility in these systems through measurements of the spin lattice relaxation rate, 
T^'^, and the Knight shift, K, of the In, Ga or M nuclei. Since the nuclear spins are 
coupled to the electron spin moments in the system, they are sensitive to the dynamics 
of these electron spins. In fact, T^^, is sensitive to the imaginary part of the dynamical 
electron spin susceptibility, whereas K probes the static susceptibility. The details of 
how the nuclei couples to the electron spins, however, are important for interpreting the 
Tf ^ and K data. For example, if a nucleus is coupled to more than one electron spin, 
then the hyperfine field can vanish for particular spin configurations at certain nuclear 
sites, depending on the symmetry of the site. 

In this manuscript we attempt to provide a complete analysis of the hyperfine 
couplings to the various NMR-active sites in this 115's. Since their discovery six years 
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Figure 1. The structure of the CeMIus materials. The Ce (Wykyoff position la) 
atoms are yellow, the In(l) (Ic) atoms are orange, the In(2) (4i) are green, and the M 
atoms are blue (lb). Details of the structural parameters can be found in |l()j . 

ago, several groups have reported Knight shift and spin lattice relaxation rate data in 
these materials. However, in order to make any connection between the NMR data and 
physically relevant quantities such as the dynamic or static susceptibility, one needs to 
understand the hyperfine interaction. Shortly after the first NMR experiments in the 
cuprates, Mila and Rice showed that the response of the Cu and O can be understood in a 
relatively simple single spin component model by elucidating the hyperfine interactions 
|12j . Their pioneering work led to significant insight into the cuprates as observed 
by NMR. The hyperfine scenario in the 115's, and the heavy fermion compounds in 
general, is more complex than that in the cuprates. There is strong evidence for two 
components of spin susceptibility, with different temperature dependences. Since the 
spin lattice relaxation is dominated by fiuctuating hyperfine fields, these two components 
can contribute in unexpected ways to the temperature dependence. Consequently, one 
must approach any interpretation of Tj~^ data with an eye to these anomalous effects, 
and a complete analysis requires an intimate understanding of the hyperfine interactions. 
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2. Transferred Hyperfine Couplings 

Most of our knowledge of the hyperfine interactions in the 115's comes from 
measurements of K and the static hyperfine field, -ffhyp in the magnetically ordered 
state. In general, a nuclear spin I is coupled to an electron spin S via an interaction. A: 

Tihyp = ihAI ■ S, (1) 

where 7 is the gyromagnetic ratio. In an external field Hq, the electron spins are 
polarized via their uniform (i.e. q=0) susceptibility: M = gfiB{S) = xHq- This leads 
to a shift, K, in the NMR resonance frequency: 

n = 'yhi{l + K)-Uo, (2) 

where K = Ax/g^is- If the susceptibility is strongly temperature dependent, then 
the Knight shift should be as well. By measuring K{T) and x{T) independently, and 
plotting K versus x with T as an implicit parameter, one can extract the value of the 
hyperfine coupling constant, A. 

In the CeMlus materials, x is strongly temperature dependent with Curie- Weiss 
behavior. There are two sets of spins that potentially can give rise to this susceptibility: 
the conduction electron spins, Sc = 1/2, and the local moments of the Ce^^ 4f^ electrons, 
Sf. Note that for the lanthanides, the spin-orbit interaction is several eV, so for 
experimentally realizable temperatures, the important quantum number is the total 
spin J = L + S, where L = 3 and S* = 1/2 for Ce^+. The degeneracy of the ground 
state 5/2 multiplet is further lifted by the crystal field interaction, which has tetragonal 
symmetry in the 115's, typically on the order of 50-lOOK, so that at low temperatures the 
ground state of the Ce^"*" is a pseudospin doublet with an effective g value determined by 
the crystal field. For simplicity, we denote the 4f spin by Sf. The local moments and the 
conduction electrons also experience a Kondo interaction and an exchange interaction 
between different 4f spins that is driven by an RKKY or other similar mechanism. 
The hyperfine interaction couples the nuclei to these two sets of spins, but the coupling 
strength is several orders of magnitude smaller than the other interactions in the system. 
Therefore, the nuclei are sensitive probes of the magnetic behavior of the system, but 
they do not significantly perturb the intrinsic phenomena. Empirically, we find that 
for sufficiently large temperatures T > T*, where T* is material dependent but on the 
order of 50K, K is linearly proportional to x- Since the dominant contribution to the 
susceptibility is from the local moments, Sf, at these temperatures, we conclude that 
the largest hyperfine field arises from the local moments. In general, however, there 
can be also an on-site hyperfine interaction to conduction electrons. A more complete 
description of the hyperfine interactions is given by: 

7ihyp = 7^I- ( A-Se+5^1-S^j , (3) 

where A is an on-site hyperfine tensor interaction to the conduction electron spin, and B 
is a transferred hyperfine tensor to the Ce 4f spins ^T] . A similar scenario is present in 



Hyperfine Couplings in CeMIn^ 



5 



the high temperature superconductors jT2]- For the CeMIus materials, the mechanism 
of the on-site couphng A is probably due to a combination of core polarization, 
and polarization of unfilled orbitals with s and p symmetry. The mechanism of 
the transferred interaction, B is probably due to a combination of orbital overlap 
between the 4f wavefunction and ligand nuclei core wavef unctions, and a polarization 
of the conduction electron spins via the Kondo interaction between Sc and Sf. For 
the remainder of this article, we assume that the hyperfine parameters are material 
dependent constants, and will not discuss their microscopic mechanism further. 



2.1. Coupling to In(l) 



In the antiferromagnetic compound CeRhlns (T/v = 3.8K) the magnetic order is an 
incommensurate spiral with the moments in the ah plane, and the modulation along 
the c direction, with Qaf = (0.5,0.5,0.297) fSJE]- In this system, the In(l) lies in 
the plane with four nearest neighbor Ce spins, and by symmetry one might expect the 
transferred hyperfine field would vanish in the ordered state. In fact, the hyperfine field 
is finite (-ffhyp ~ l.SkOe), and is about five times larger than one would expect from 
a direct dipolar interaction between localized 4f spins and the In (/ = 9/2) nuclear 
spin. The reason the transferred hyperfine does not vanish is that the eigenvectors of 
the tensor B do not coincide with the unit cell coordinates, but rather point along the 
bond directions between the Ce and the In(l). We can write B as: 
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(4) 



in the bond-coordinate system, where Ba corresponds to the transferred hyperfine 
coupling along the a direction as shown in Fig. |21 In the unit cell (and experimental) 
reference frame, this tensor becomes 
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for the sites i e (1, 2, 3, 4) in Fig. d where Bq = {B\\ + 5^) /2 and Ba = {B\\ - Bj_) /2. 

is probably anisotropic, but diagonal in the unit cell 
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The on-site interaction, j 
coordinates: 
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Figure 2. The principle axes of tlie transferred hyperfine tensors B.^ of tlie In(l) and 
the In(2) sites. 



The Knight shift at the In(l) site, then should be. 

K^{T) = + AB^{S])) /Ho (8) 

where the approximation should hold for temperatures T > T*, where the Curie- Weiss 
susceptibility of the 4f electrons dominates the Pauli susceptibility of the conduction 
electrons. Below, in section 3, we discuss the Knight shift for T < T*, where the 
approximation no longer holds. 

In the antiferromagnetic state of CeRhlus the ordering is given by: 

(S/(r)) = 5*0 cos ( — j cos ( — j (cos(27rgo-2)x + sin(27rgo-2)y) , (9) 

where go is the wavevector of the spiral along the c direction jj^l El- The hyperfine 
field at the In(l) site is then given by: 

Hhyp = 45aS'o (sin(27rgo2;)x + cos(27rgo2;)y) • (10) 

ifhyp oc — would vanish if B did not have off-diagonal terms in the coordinates 
of the unit cell. Also, since Hhyp -L c, and the electric field gradient (EFG) vector 
for the In(l) is parallel to c, the resonance frequency in the antiferromagnetic state is 
independent of z. Consequently, even though the magnetism is incommensurate, the 
In(l) lines remain sharp [13j. The orientation of the 4f moments and the hyperfine fields 
are shown in Fig. El 

A further consequence of anisotropic transferred hyperfine tensor, B is that the 
ln(l) is sensitive to the critical fluctuations of the order parameter above T/v- We will 
come back to this subject later in a discussion of the form factors. 
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Figure 3. The hyperfine fields in the CeMIns materials. The atoms are the same as 
in Fig. n The arrows on the Ce atoms denote the direction of the electronic moments, 
for the ordered structure in the CeRhlns antiferromagnetic state. The arrows on the 
other atoms (In(l), In(2) or M) denote the direction of the hyperfine fields at the 
particular site, as discussed in the text. An animated version of this figure is available 
online, showing how the hyperfine fields change as the direction of the Ce moments 
changes. 



2.2. Coupling to In(2) 

The In(2) is in a low symmetry site out of the Ce plane. We expect that the transferred 
hyperfine coupling tensor to this site should also be diagonal along the bond axes (see 
Fig. 121). In the unit cell coordinate system, we have then: 

/ Bo + B„cos2e 5„sin2^ \ 
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Ba COS 20 



B^ sin 26 
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(11) 
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and 



f Bo + Bacos2e 
R, 



v 
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where 9 is the angle between the Ce-In(2) direction and the a axis. For CeRhlns 9 
In this case, the Knight shift is given by: 



Kb 



2{Bo + BaCos{29)){S'})) IH, 
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(12) 

45°. 

(13) 
(14) 
(15) 



A,{S:) + 2{Bo - BaCos{29)){S})) /Ho, 

where we have ignored the small component of hyperfine field that is perpendicular to 
Ho- As in the case for In(l), we assume that the Pauh susceptibility is negligible for 

T > T\ so 



Ka 



(16) 



K, 

Kr 



2{Bo + BaCos{29))xab{T) 
2ByXab{T) 

2{Bo-Ba cos{29))xc{T). 

Even though the spin susceptibility is isotropic in the ah plane, the In(2) shift depends 
on the whether the field is in the a or 6 directions. The reason for this is that the 
symmetry of the In(2) site is lowered when the field is in the plane, so two of the In(2) 
sites have the field applied parallel to the face of the unit cell, and two of the sites have 
a field perpendicular to the face of the unit cell. 

In the Neel state, the hyperfine field at the In(2) site is parallel to the c axis: 

Hhyp = 25aS'oCos(27rgo2;)z. (17) 

Once again, the hyperfine field at the In(2) site is finite only because the tensor B is 
not diagonal in the unit cell coordinates, and thus the In(2) is also sensitive to critical 
fiuctuations. The orientation of the hyperfine field at the In(2) sites is shown in Fig. 
El Note that in the case of the In(2), the resonance frequency is dependent on the 
z coordinate, therefore the In(2) exhibits a broad "powder pattern"-type resonance 
spectrum in the ordered state of CeRhlns [T3j. 

2. 3. Coupling to M 

For the M site, the bond axes coincide with the unit cell axes, so the hyperfine tensor 
is diagonal in the unit cell basis. We thus have: 

/ Bab ^ 

(18) 
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and 



Kat 



'MS:) + B,{S})')/Ho, 



(19) 



(20) 



so for T > T* 



Kab 



2BabXab{T) 



(21) 



In the interest of notational simplicity, we have not included any indices on the on- 
site and transferred hyperfine coupling constants, and B^ for the different sites 
(In(l), In(2) and M). In general, each one of these couplings is different, and is material 
dependent. 

3. The Knight Shift Anomaly 

In a generic Kondo lattice system, there are two spin species, Sc and Sf. Therefore, there 
are three distinct spin susceptibilities: Xcc = (ScSc), Xcf = (ScSf), and Xff = (SfSf). 
In Eqs. Uni and we made the approximation that {Sc) ^ (5'/) for T ^ T*. A 
more complete description is given by gcfJ^siSc) = {Xcc(T) + Xcf{T))Ho and gffisiSf) = 
{Xff(T) + Xcf (T)) Ho, where gcj are the g-factors. In fact, the full expression for the 
Knight shift is given by: 



where rj refers the the particular nuclear site, a is the field direction, and for simplicity 
we have absorbed the g-factors into the definition of the hyperfine constants JI] . The 
bulk susceptibility is given by: 



Note that if Aa = X^ienn ' then the relation i^' oc x is recovered, however in general 
this relation will not hold. If Xcc(T), Xcc(T), and Xcc(T) have different temperature 
dependences, then the Knight shift will not be proportional to susceptibility, leading to 
a Knight shift anomaly at a temperature T*. 

Since a complete solution of the Kondo lattice remains a challenge for theory at 
present, there is no description of the temperature dependence of Xcc(T), Xcc(T), and 
Xcc{T)- However, using empirical observations Nakatsuji and coworkers developed a 
phenomenological two-fluid picture of the susceptibility and specific heat in CeCoIns 
which also works well for understanding the Knight shift anomaly fSlElini- In this 



picture, Xcf{T) ~ (1 - T/T*) \og{T*/T) for T < T*, and Xcf{T) ~ for T > T*, and we 



assume Xcc is negligible for all temperatures. This relation appears to hold for several 




(22) 
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(23) 
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different heavy fermion and mixed valent compounds down to the relevant ordering 
temperatures [TT[[Tn]. Recently we have found evidence that XcfiT) saturates below a 
temperature Tq -C T* in CeCoIns for H > Hc2- Recent theoretical work lends support 
for the two-fluid picture [TTl HH] , but a full description is still lacking. Nevertheless, the 
NMR data in a wide range of materials strongly suggest that a complete theoretical 
description should involve one localized degree of freedom giving rise to Sf, and another 
delocalized degree of freedom, giving rise to Sc- 



4. The Spin Lattice Relaxation Rate and the Form Factors 

The spin lattice relaxation of the nuclei in the CeMIns materials is dominated by 
magnetic fluctuations mediated by the dynamic hyperfine fields at the nuclear sites. 
The spin lattice relaxation rate, T{"^, can be written in general as: 

1 r°° 

7f^ = ^nJ2 {mwo)y^''dt, (24) 

J-l,a p Jo 

where we have dropped the "hyp" notation for simplicity, ujl is the Larmor frequency, 
jn is the nuclear gyromagnetic ratio, and the sum is over the two spatial directions 
orthogonal to the a direction. In general, the hyperfine field is given by: 

H(t) = A ■ S"(0, t) + J2 ■ S^(ri, t) (25) 

ignn 

where the i-sum is over the nearest neighbor Ce spins, and the 7 sum is over the spatial 
directions. This expression can be simplified by writing: 

Hp{t) = 5^F|^(q)5^(q,t) + F/^(q)5^ (q, t), (26) 

q>7 

where Sj^''^\ci,t) = "^j. Sj^''^\r,t)e~'^'^''^ , and the form factors F('^'''^)(q) are given by: 

F|^(q) = Ap^6p, (27) 

^/7(q)= (28) 

i 

We then have: 

{H,{t)H,m = 5^(F|^(q))^SW(q,t)4^)(q,0)) + 



q.7 



2 E ^l7(q)^/7(q)(4^nq, t)s\^\^, 0)) + 

q,7 

E(^/,(q))'(5ynq,i)5ynq,o)). (29) 

q.7 

Here we have assumed that spin correlations at different q— vectors and directions are 

uncorrelated: {Sj^''^\q^, t)Sy''^\c{', 0)) = Sq^q;6f^^{Sj^''^\q^, t)Sj^''^\q^, 0)), although such an 

assumption may not be always valid [12] . Using the fluctuation-dissipation theorem, Eq. 

121] can be rewritten as: 

1 _ 1 2 1 
^ - ^ + + — T' (^°) 

JiQ 
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where: 



,„,,-7?.*.riim$:*r(q)i^. (31) 

l,a q,/3 

where x^'r?'/?!^' ^) the dynamical spin susceptibihty along the /5 direction, and 

0r'(q) = E^/5?(q)^/57'^(q)- (32) 

7 

Clearly, the behavior of the spin lattice relaxation predicted by Eq. 1201 can be 
quite complex, depending on the temperature dependence of the various dynamical 
spin susceptibilities Xcc x'cf X//, especially if the temperature dependence of these 
three quantities differ. At present, there is no clear theoretical model for how these 
susceptibilities should behave. However, we can make several important and general 
conclusions. NMR Knight shift studies and theoretical considerations suggest that 
Xcf is negligible for T > T* [201 ■ Therefore, we expect that Tf^ ^ T{'L + T^jf for 
T > T* . Below T*, a new contribution, T^l^ can become important, and may lead to 
an anomalous temperature dependence, as observed in CeCoIn5 for fields applied in the 
ah plane pTj . 

The filter functions (q) given in Eq. ESI enhance or suppress the relative 
importance of spin fluctuations to relax the nuclei at different parts of the Brillouin 
zone. For the In(l) we have: 

0f^(q) = lQBlcos\q^a/2)cos\qya/2) 

+ 165^ sin2(g:,.a/2) sm^{qya/2) (33) 
0//(q) = 165^2 cos2(g^a/2)cos^(gya/2) (34) 

Note that -F|(q) = Aj^p is q— independent. For the In(2), 

0{^(q) = 4(5o + 5„ cos 2^)2 cos2(g,a/2) + 

4(5„ sin 2^)2 sin2(g^a/2) (35) 
0f(q)=45jcos2(g,a/2) (36) 
0{^(q) = 4(5o - Ba cos 2ef cos\q,a/2) + 

4:{BaSm2efsm'^{q^a/2) (37) 

and for the M site: 

0{^(q)=45^,cos2(g,c/2) (38) 
0f(q) =4S2cos2(g,c/2) (39) 

These quantities are plotted in Figs. lU and El The anisotropic nature of the hyperfine 
tensor of the In(l) and In(2) gives rise to a second component in their filter functions 
that does not vanish at q = (tt/o, ir/a, q^). A priori, one might expect that the In(l) and 
In(2) would be insensitive to fluctuations at these wavevectors because they are located 
in symmetric positions in the unit cell with respect to the Ce. However, since the bond 
axes are not coincident with the unit cell axes, the hyperfine tensor is not diagonal. As 
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Figure 4. The filter function <pjf,{(l) for the In(l), shown for Ba/Bo — 0.6. Note that 
(/>{{(q) remains finite for q = (tt/o, n/a, qz). 



a result, the In(l) and In(2) are sensitive to antiferromagnetic fluctuations ^2T]. In the 
cuprates, there is no such anisotropy in the hyperfine tensors, since the bond axes of 
the Cu and the O nuclei lie along the unit cell directions. 

It should be pointed out that the the spin lattice relaxation at the In(2) site is 
particularly difficult to interpret, because the electric field gradient (EFG) does not 
have axial symmetry. The quadrupolar part of the hamiltonian of the In(2) is given by: 



rio 



6 



(40) 



where the the eigenvalues of the EFG tensor j22]. Since Tig does not commute 

with la, h or /c, the eigenstates of the nuclear levels are non-trivial superpositions of 
the |m) states, and the usual selection rules for transitions between nuclear levels do 
not hold. Consequently, unlike the In(l) which does have axial symmetry, the master 
equation which governs the relaxation of the In(2) nuclei cannot be solved analytically. 
One can measure the relaxation at the In(2) site, which is a function of the rates given 
in Eq. but it is not straightforward to independently measure the rates T^^l. 
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Figure 5. The filter functions 4'ff{ci) and (f)l'^ (q) for the In(2), shown for Ba/ Bq = 0.3. 
Note that (f>{/ {(l) remains finite for all q, whereas ipl^ {(\) does not. 

In fact, several groups have measured the temperature dependence of the relaxation 
rate at the In(l) site in the various 115 compounds, and have found a rich variety 
of behavior both in the superconducting and normal states [23 123 1211 I2S1 I2Z|- In 
CeRhlns, the temperature dependence of the spin lattice relaxation is a strong function 
of pressure j2ni , and even shows signs reminiscent of the pseudogap observed in the high 
temperature superconductors. One explanation is that the spin fluctuations are modified 
with pressure as the ground state evolves from antiferromagnetic to superconducting. 
However, the internal field in the ordered state is reduced linearly with pressure, while 
the ordered moment remains constant [2H]- This result suggests that the hyperfine 
constants may change with pressure, a result that has not yet been confirmed by Knight 
shift measurements under pressure. If the hyperfine constants are pressure dependent, 
then the temperature dependence of the dominantly reflects these changes, rather 
than changes in the spin fluctuation spectrum. Measurements of T^^ in CeRhi_a;Ira;In5 
alloys and in CeRhlns under pressure suggest that there are regions of the phase 
diagram where both antiferromagnetism and superconductivity coexist microscopically 
|23| l27] . In such cases, it appears that the typical temperature dependence in the 
superconducting state of d-wave superconductors may be modified, possibly due to the 
presence of the antiferromagnetic order. Further studies of the Knight shift in these 
systems may shed important new light on this behavior. 

5. NMR in the Actinide 115's 

Recently the actinide based 115's, synthesized with U, Np and Pu have attracted 
attention due to the large superconducting transition temperature in PuCoGas [HE]. 
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One of the striking discoveries in PuCoGas is that the temperature dependence of 
spin lattice relaxation rate scales with the transition temperature in the same manner 
in CeCoIns, PuCoGas, and YBa2Cu307 suggesting a common mechanism. Although 
there is little detailed information available about the hyperfine couplings present in the 
PuMGas materials, there is no reason to suspect that the general form of the tensors 
A and B should differ from those in the CeMIus materials. Therefore, the spin lattice 
relaxation rate measurements, whether on the In(l), In(2), Ga(l) or Ga(2), should be 
sensitive to antiferromagnetic fluctuations at (vr/a, vr/a, g^). The unusual temperature 
dependence of T^^ may well reflect the temperature dependence of X// and Xcj- The spin 
lattice relaxation rate measured in PuRhGas, however, differs dramatically from that 
found in PuCoGas below ~30K, well above T^. j^H]- Although this anomaly may reflect 
unusual spin dynamics jHO], another explanation may be the presence of a Knight shift 
anomaly at 30K observed in this material [SI] . This discrepancy highlights the unusual 
hyperfine interactions present in all of the heavy fermion materials, and demonstrates 
the need for a general understanding of the two spin components and their dynamics 
before placing too much emphasis on the behavior of the spin lattice relaxation as a 
function of temperature. 
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